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Vy, 
ABSTRACT 


Areas of convective precipitation were delineated 
and averaged convective precipitation amounts determined 
within these areas. A statistical, step-wise, screening, 
linear regression procedure was used to correlate certain 
large scale parameters with.the convective precipitation 
amounts, 

Of the parameters considered in the two cases 
investigated, the best predictor of convective precipita- 
tion was either the terrain-induced vertical motion at 
the top of the boundary layer, or the product of the 


latter and the specific humidity at the same level- 
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1. INTRODUCTION 
The failure of numerical models to account for 
convective precipitation is well known (see, e.g., 
Danard (1966 a, 1966 b)). Such a shortcoming is expected 
Since the equations used for numerical prediction are not 
applicable to sub-synoptic scale motion. Convective 
precipitation is highly variable in both time and space. 
This would prevent such precipitation from being accu- 
rately accounted for even with a much denser observational 
network than the one available in the United States today. 
This investigation is an attempt to delineate the 
areas of convective precipitation by the use of semi- 
empirical methods. Convective precipitation will be 
interpreted here as the average amount associated with 
Ssub-synoptic scale systems over a Square area whose side 
is equal to the distance between adjacent grid points 
used in numerical prediction (381 kilometers at 60 degrees 
north). The amount defined is assumed to be determined 
by the large-scale properties of the atmosphere. 
Convective activity is confined to areas of unstable, 
or conditionally unstable, lapse rates. That is, the lapse 
rate should be steeper than the moist adiabatic lapse 
rate. Kuo (1965) describes a convective cloud model which 
will "bring surface air to all levels up to a great 
height so that inside the cloud the vertical distribution 


of temperature and mixing ratio are those of the moist 


adiabat through the appropriate condensation level.” 

This suggests that the rate of production of cloud and 
total water content would depend on the specific humidity, 
or mixing ratio, and the vertical velocity at some lower 
boundary. Ooyama (1963) assumed that the rate of 
precipitation was proportional to the flux of water 

vapor through the top of the frictional boundary layer. 
The total gain of water vapor in an air column, assuming 


no flux through the top of the column, may be represented 


by 


P. 
T= GleVe Wap — (wnqn)/g 5 9 (1) 


(e.g., Krishnamurti (1968)), where I is the water gained, 
q is the specific humidity, and p is pressure. The sub- 
scripts.b and t indicate the bottom and top of the column, 
respectively. Here the first term represents the flux of 
water vapor through the sides of a column of unit area, 
and the second term represents the flux through its base, 
which 18 taken here as the top of the planetary boundary 
layer. Since specific humidity decreases rapidly with 
height, the second term is frequently larger than the 


first if the level Py is suffietbently hiteh. 
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2. PROCEDURE 

For the two case studies included here, Danard 
(1966b) has made numerical predictions of gridpoint 
precipitation within an area shown in Figures 5-8. 
Predicted amounts obtained from the first 12 hours of 
these integrations were assumed to represent precipitation 
effects of large-scale motions. Observed precipitation 
amounts were obtained by averaging all individual values 
reported by synoptic stations over a Square centered 
on the grid point and having a side equal to the distance 
between adjacent grid points (approximately 381 kilometers). 
Convective precipitation amounts used are the difference 
between the averaged observed amounts and the amounts 
predicted by Danard. Only those points where the lapse 
rate was close to, or greater than, the pseudo-adiabatic 
value (see below) were used. The maximum 12-hour 
precipitation amount predicted by the numerical scheme 
in these areas was only 0.4 millimeters in either case 
study. Thus, the areas of convective precipitation were 
clearly delineated from those of large-scale precipitation. 

A statistical, step-wise, linear regression procedure 
was used to correlate convective precipitation with the 
following parameters: (1) the difference between the pseudo- 
adiabatic and the computed lapse rates (R,-[), (2) the 
vertical velocity (Ww) at the top of the boundary layer, 


(3) specific humidity (q) at the top of the boundary 
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layer, (4) the product of (2) and (3), and (5) dew point 
depression (Tels) in the lower troposphere. The program 
used is available in the Naval Postgraduate School 
computer library, where it is listed as BMDO2R. A 
sequence of multiple linear regression equations is 
computed in a step-wise manner. At each step the 
variable is added which makes the greatest reduction 

in the error sum of squares. 

For the two cases considered, average temperatures 
and dew-point temperatures were computed for the 850, 
700, and 500 mb levels. The reported values at the 
beginning and end of the period were used to compute 
time averages. These averaged temperatures were used 
in the computations. 

All averaged parameters were computed from values 
reported by the stations. The resulting fields were then 
analyzed, and the values for the grid points used were 
taken from the analysis. 

As used here, the lapse rate is defined to be the 
difference between the temperatures at 850 and 500 mb. 
The pseudo-adiabatic lapse rate between 850 and 500 mb 
was taken from an Arowagram as a function of the tempera- 
ture at 700 mb. Convective precipitation was considered 
possible where the difference between the lapse rates 
- (M.-P) was less than 2.0°C/(350 mb). The areas delineated 
by this criterion included all regions where large amounts 


ef convective precipitation was observed in both cases. 


dee 


The dew point depression used was averaged between 
850 and 700 mb. Here the 500 mb level was omitted 
because, in both cases many stations reported motor- 
boating at that level. In the cases where motor- 
boating was reported at the lower levels, the minimum 
dew point depression for which a radiosonde normally 
reports motor-boating was used. 

Danard (1966a) expressed the vertical velocity at 


the top of the boundary layer as 


W = Wr + We 5 (2) 


where Wy and we, are the vertical velocity induced by 
orographic and eddy-stress influences, respectively. 

In the area considered in this study (see Figure 1) 
most of the terrain heights are less than 400 meters 
above sea level. Therefore, the surface geostrophic 
wind was taken to be that at 1000 mb. Thus, Wy is well 


approximated by 


Wr = (Vz ooo # VP ? 37) 
where Py is the standard atmosphere pressure at the 
elevation of the smoothed terrain provided by Berkofsky 


and Bertoni (1955), and ( V_) 


¢/1000 is the geostrophic 


wind at 1O0OO mb. 
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The frictionally induced vertical velocity was 


computed by 
Wp = IAF(KVXCIVYGVY)  . (4) 


Here V. is a representative geostrophic wind in the 
planetary boundary layer and Cy is the drag coefficient 
(Cressman (1960)). In the cases considered here ( V_) 
ge“ L000 

was used. However, over higher terrain, or in cases 
where ( V.) 1000 cannot be easily computed, ( V,.) 850 
might be more useful. 

Specific humidity at the top of the boundary layer 


may be expressed as 


= 620 
q = peti ) (5) 


where e is the vapor pressure and p is the pressure in 
mb at the top of the boundary layer, assumed to be 850 
mb here. Vapor pressures used were the saturation 
values (Smithsonian Institution (1951)) corresponding 


to the dew point temperatures at 850 mb. 
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3. RESULTS OF CASE STUDIES 

Two synoptic situations were investigated: 1200 GMT 
11 February to 0000 GMT 12 February, 1965, and 1200 GMT 
16 March to 0000 GMT 17 March, 1965. These intervals will 
be referred to as Cases A and B, respectively. The 
surface synoptic situation for the beginning and end of 
the two 12 hour periods may be seen in Figures 1-4. 

Both cases were characterized by deepening cyclones 
moving from the southwest. These storms produced some 
of the worst weather of the 1965 winter in the central 
United States. 

The number of grid points used in Cases A and B 
(i.e., the sample sizes) was 16 and 17, respectively. 
Tables 1 and 2 give the simple correlation coefficients 
between all the variables used for Cases A and B, 
respectively. In both cases the correlation between w 
and Wq was very high. In Case A, wq had the highest 
simple correlation with convective precipitation (P). 

For Case By the correlation between wand P.was slightly 
higher than that betweenwq and P. The second highest 
simple correlation for Case A was between (T-T, ) and P. 
However, there was little or no correlation between these 
variables in Case B. 

The results of the step-wise regression analysis may 
be seen in Tables 3 and 4 for Cases A and B, respectively. 


In the step-wise regression procedure, the first variable 
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entered is always the one with the highest simple 
correlation with P, unless the "deletion option” of 
the step-wise regression program is invoked. In Case 
A, after wq entered the regression equation, w was 
essentially redundant and did not enter the equation 
at all. For the same reason wq did not enter the 
equation in Case B, since W entered first. In both 
cases there was little improvement in the multiple 
correlation between the predicted value and P after 
the first variable entered the regression equation. 

An attempt was made to determine if a time delay 
between the flux of moisture through the boundary 
layer and resulting precipitation could be detected, 
To accomplish this, the regression procedure was 
repeated using variables computed for the beginning, 
rather than for the average, of the time period. The 
resulting»simple correlation coefficients are given Pn 
Table 6 and 7. In Case A, there was noticeable improve- 
ment in the correlation between both wq and P, and Ww 
and P. However for Case B there was little change in 
these corre lari Ons . 

The step-wise regression program used permits the 
variable wq to enter the regression equation first for 
Case B if the variable Wis deleted. In this way it is 
. possible to obtain the variable wq as the primary 


predictor for both Cases A and B. Ina similar manner yw 


16 


has been selected as the primary predictor for Case A. 
The resulting equations, with only one variable entered, 
may be seen in Table 5. Due to the high simple corre- 
lation between W and wq there is little loss in explained 
variance for either case, when W or wq is employed as 
the primary predictor. 

In order to improve the prediction equations and 
to ensure that P2O, predictions were made for both 


cases using the equation 


ul 


Ce Ae + A,(q) , for wa <-Ad/A, (6) 


=) 5 for wg Z—-Ad/A,y a 


The coefficients &. and A. were determined by computing 


I 
the quanity 


S?=Z(P,-PY (7) 
where the summation extends over all the points in the 
Sample. P 5b is the observed convective precipitation at 
a gridpoint and P is the amount given by Equation (6). 
se is Simply the error sum of squares. Various values 


Ole A, and A, were tried until the minimum in ge was 


1 
found. This procedure was repeated using w& as the 
predictor, The improved prediction equations are shown 

in Table 9. Comparing the last columns of Tables 5 and 

9, it is seen that the above procedure results ina slight 


reduction in error for Case B. However, for Case A the 


improvement is negligible. 
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Predictions were made for both cases using the 


equation 


P = -0.45 - 1.20(wWq) , for wa e- 0.375 


(8) 
= © Sfor Wak-OU57> . 


Equation (8) was obtained by averaging the coefficients 
given"in the first two lines’ of Table @®. The Yresuilies 


may be compared with the observed convective precipitation 


in Figures 5-8. 
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TABLE I 


Matrix showing the simple correlation 
coefficients between the variables for 
Case A, P is convective precipitation. 











E oa i uw q wq T-T5 

Pp I. O@ 0.28 -0.43 O.l2 “0.51 -O.47 

wot 1.00 -O, 34 -0.13 -0.35 -0.44 

UW ie 00 -0.10 0.95 O55 

q 1.08 -0.28 OmOl 

ae I GO Om 3 
TAB EEL 1 


same aS Table I for Case B 


EP Pie I w) q wd T-T 





P i OG -0.25 -0.67 -0,24 -0.65 = 
noe 1.00 0.01 -0.65 -O.01 -0.34 
u) 1.00 0.45 . 0.98 O531 
q 1.00 0.48 0.36 
wd 1.00 0.28 





TABLE III 


Results of step-wise regression analysis for Case A in the 
form of regression equations. p= AS Se ACV. , where Peis 
the precipitation, in mm/12 hr, and V, is the variable 
entered at the i step. R is the multiple correlation 
coefficient between P and the Ppredictand at the i step. 


| «1. [SPI Tia, .,,,— 


Root mean 
Square error 


Variable Regression coefficients of estimate 
Step entered Wits R a ag A, A, Ay (mm/12 hr) 
= as 
il wWwq dynes em “sec O ea 1 Os ay -1.66 4.65 
2 waa SC 0.56 4.73 Eilgege See 4 48 
3 R7P  °c/350mb 0.56 4.47 -1.17 -0.39 -0.30 447 
d, q em/keg O56 3.73 or. 14 -0,39 mO.35 -O.08 64.47 


NOTE: wq is multiplied by 10° 


tt °®©|+O &2« ©. OOO 
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TABLE, IV 


Same as Table tpl tor Cases: 


Root mean 


square ‘error 
Variable Regression CoetMicienys of estimate 
entered Uns ts R A Ay A, Ax Ay (mm/12 hr) 


















Step 
i w dynes geo” 0.67 -0.08 S15 27 I te, 
2 mor °c /350mb 0.72 -0.33 -1.36 -0.17 age. 
3 q em/kg e. 7 wees mi 12 ee Once) 0.52 mW, 1e 
mn TT zc Ona BO. 70 Mae «-Ouee. -=0S1 «(O04 las 
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TABLE V 


Regression equations using only 
tne first variable entered. 

















Regression Root mean 
Variable coefficients Square <emror 
Case entered R A. A, of estimate 
A wd Onl Oy 57 -1.66 4.65 
B wg 0.65 -0.01 -0.32 i 
A w 0.43 LTyO7 =e 4.89 
B ud) 0.67 =@F@s.  -1i37 1.28 
TABLE VI 
Matrix showing the ssimele correllavion 
coefficients between the variables for 
Case A. The primed variables were computed 
for the beginning of @the time period. 
P n-r uw)’ q/ wa! T-T) 
iP 00 0,23 -0.65 0.04 -0.64 ~ORbO 
Tne 1.00 -0,24 -0.39 -0.14 -OF50 
uw" 1.00 -0.40 0.99 0.36 
q’ 1 ee -O,42 -0.62 
wg! 1.00 O43 





TABLE VII 


Same as Table VI for Case B 


| me 


[ wn ud)’ q! w'q’ Te, 

P 1.00 -0.25. 0.62 -O0. 41 -0.66 @e37 

vig 1 260 0.06 -O.644 0.64  -0.o 

vu’ 1.00 0, 44 @nos) Geese 

qr 1.66 0.41 0.01 

wa! 1.00 @.26 
TABLE VIII 


same as Table V for variables 
computed fer the begiantug wr 
the time period. 





Regression Root mean 
Variable coefficients Square error 
Case entered R a Ay of estimate 
A wa!’ 0.64 1.71 #1£x°-2.08 Lae 
B a)’ a! 0.66 O821 70.728 1 Sh 
A ww’ 0.65 2.22 -13.77 44 : 
B u)’ 0.62 6.07 =2.59 1.36 





ae) 


Case 


a = WwW & 


Improved regression equations using 
only the first variable entered. 


Variable 
entered 


Ww q 
wg 
uw) 
w 


TABLE IX 


Regression 
coefficients 
i Ay 
0.20 =1 80 
-~1.10 -O0,50 
Voo -0.35 
-0.80 -1.85 
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Root mean 
Square error 
of estimate 


4.63 
1 ly 
4,88 
Tree 
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Figure § 


Computed 
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Figure 6 
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Figure 7 
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Precip. 





Figure 8 
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4, CONCLUSIONS AND REMARKS 

Probably the most that could be claimed for the 
results of this study is that it may have some 
applicability over the central and eastern United States 
for the winter months during the period 1200 to OOOO GMT. 
However there is indication that the product of vertical 
velocity and specific humidity at the top of the 
boundary layer can be used to predict convective 
precipitation. 

While this study is diagnostic in nature, it uses 
data that can be predicted numerically. Consequently, 
the technique could be included in a™prognostic scheme 
with some expectation of success. The technique of 
using lapse rates, as described here, to determine areas 
of possible convective activity could be applied to 
local forecast areas using forecast temperature fields. 

For a synoptic situation to be used in this type 
of a study, there must be a clear delineation between 
large-scale and convective precipitation. To find out 
if a situation meets this requirement, a prediction of 
precipitation amounts due to large-scale processes must 
be made and compared to the reported precipitation. The 
time required for this procedure limited this study to 
the two cases considered here. Furthermore, the small 


Size of the samples virtually precludes a statistical 


es 


proof of the significance of the results. However, 
the fiindimgs-ofethis study are amffieienthy encoun@ging 


to justify ferthem work. 
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